Future constraints on variations of the fine structure constant from combined CMB 

and weak lensing measurements 
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We forecast the ability of future CMB and galaxy lensing surveys to constrain variations of the 
fine structure constant. We found that lensing data, as those expected from satellite experiments as 
Euclid could improve the constraint from future CMB experiments leading to a Act/ a — 8 x 1CP 4 
accuracy. A variation of the fine structure constant a is strongly degenerate with the Hubble 
constant Ho and with inflationary parameters as the scalar spectral index n s . These degeneracies 
may cause significant biases in the determination of cosmological parameters if a variation in a as 
large as ~ 0.5% is present at the epoch of recombination. 



I. INTRODUCTION 

After the recent measurements of Cosmic Microwave 
Background (CMB hereafter) anisotropies, galaxy clus- 
tering and supernovae type la luminosity distances (see 
e.g. [l|-y|) cosmology has entered a "golden age" where 
not only a " standard" model of structure formation has 
been firmly established but also where physical assump- 
tions not testable in laboratory can now be probed with 
unprecedented accuracy. 

One of the most promising aspects of this new cos- 
mological framework is indeed the possibility of testing 
variations of fundamental constants on scales, times and 
energies drastically different from those on earth. 

A variation of fundamental constants in time, while 
certainly exciting, represents a radical departure from 
standard model physics. In the past years much of 
attention has been focused on the fine structure con- 
stant a mainly because of the observational indication 
of a smaller value in the past, at cosmological redshifts 
z = 0.5 — 3.5, from quasar absorption systems data with 
Aa/a = (-0.72 ± 0.18) x 10~ 5 (Q). More recently, 
these observations have been complemented with new 
data from the Very Large Telescope (VLT), covering a 
different region of the sky. This data shows an oppo- 
site evolution for a, i.e. an increase respect to the local, 
present value, by approximately the same amount (Q) 
and leading to the hypothesis of an anisotropic variation 
in a. While all these measurements could be affected by 
some hidden systematics the search for variations of the 
fine structure constant in the early universe is clearly a 
major and timely line of investigation in modern astro- 
physics (see e.g. @ for a review). 

CMB anisotropies have provided in the past years 
a powerful method to constrain variations of the fine 
structure constant. Since the CMB anisotropies were 
generated at the epoch of recombination, approximately 
~ 340.000 years after the Big Bang, they probe the value 
of the fine structure constant at that time, when the uni- 
verse was nearly isotropic and homogeneous. Constraints 
have been obtained analyzing CMB data (see e.g. @- 



[T3|) with an accuracy at the level of ~ 10 — 1%. In the 
most recent analysis, parametrizing a variation in the fine 
structure constant as a/ao, where ao = 1/137.03599907 
is the standard (local) value and a is the value during the 
recombination, the authors of [l5[ found the constraint 
a/ao — 0.984 ± 0.005, i.e. hinting also to a more than 
two standard deviation from the current value. 

While the current CMB bound is considerably weaker 
than the observed variation from quasar spectra, the 
CMB recombination occurred at a time period corre- 
sponding to redshift z ~ 1300. It is quite possible that 
a has larger variations at higher redshifts, i.e. there is 
no reason for the variation to be constant in time. The 
CMB bound provides therefore important constraints on 
the running of the fine structure constant. 

In the next few years, we expect a further significant 
improvement in the quantity and quality of cosmological 
data. The Planck satellite mission (see [HI), for ex- 
ample, will provide accurate temperature CMB maps by 
the end of this year. On the other hand new and larger 
galaxy surveys are either already operating either under 
study. Some of these surveys will provide new galaxy 
weak lensing measurements that, when combined with 
Planck, will drastically improve the constraints on cos- 
mological parameters. The Euclid satellite mission 
selected as part of ESA Cosmic Visions programme and 
due for launch in 2019, probably represents the most ad- 
vanced weak lensing survey that could be achieved in the 
current decade. 

Future weak lensing surveys will measure photometric 
redshifts of billions of galaxies allowing the possibility a 
tomographic reconstruction of growth of structure as a 
function of time through a binning of the redshift distri- 
bution of galaxies, with a considerable improvement of 
cosmological information (e.g. dark energy plj ]; on neu- 
trinos [22l l23ll; the dark matter distribution as a function 
of redshift 26;] and the growth of structure [24j, |25[ ) . As 
far as we know, there is however still no study in the lit- 
erature that considered the gain in constraining variation 
in fundamental constant from these surveys. 

It is therefore timely to forecast the constraints on vari- 
ation of the fine structure constant achievable from weak 
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lensing data from the Euclid satellite mission. In this pa- 
per we indeed perform this kind of analysis. As expected, 
weak lensing probes are shown to be complementary to 
CMB measurements and to significantly improve the con- 
straints on variation in the fine structure constant. 

The structure of the paper is as follows. Scctionlnlcon- 
tains the description of the simulated future data used in 
our analyses. The results from our Markov Chain Monte 
Carlo (MCMC) analyses are presented in Sec. Mil We 
draw our conclusions in Sec. IIVI 



II. FUTURE DATA 

In this section we describe the method adopted to sim- 
ulate the future CMB and weak lensing surveys data used 
then to forecast the constraints on a and the remianing 
cosmological parameters. 

The fiducial cosmological model assumed in producing 
the simulated data is the best-fit model from the WMAP 
seven year CMB survey (see Ref. The parameters 

are: baryon density il^h 2 = 0.02258, cold dark mattter 
density Q c h 2 = 0.1109, spectral index n s = 0.963, optical 
depth r = 0.088, scalar amplitude A s = 2.43 x 10" 9 and 
Hubble constant Ho = 71. For the fine structure constant 
we assume either the standard value a/ao — 1, either a 
small variation a/ao — 0.996, that we believe could be 
detectable with these future data. 

As stated in the introduction we consider CMB data 
from Planck and galaxy weak lensing measurements from 
Euclid. For CMB data the main observables are the Ci 
angular power spectra for temperature, polarization and 
cross temperature-polarization. For weak lensing data 
we consider the convergence power spectra P(£) following 
the procedure described in [27| . All spectra are generated 
using a modified version of the CAMB code [28[ taking 
into account the possible variation in the fine structure 
constant as discussed in [13. 



A. CMB data 



For each channel we consider a detector noise of it; 1 = 
{9a) 2 , with 9 the FWHM (Full- Width at Half-Maximum) 
of the beam assuming a Gaussian profile and a the tem- 
perature sensitivity AT (see Tab. U for the polarization 
sensitivity). To each Ce fiducial spectra we add a noise 
spectrum given by: 

Ni = w- 1 exp(e(t+l)/%), (1) 

with 4 given by £ b = V8 In 2/9. 

Alongside temperature and polarization power spectra 
(Cf T , Cf E and Cj E ) we include also the the deflec- 
tion power spectra Cf d and Cj d obtained through the 
quadratic estimator of the lensing deflection field d pre- 
sented in HI 

d{a,b)¥=ni ]T W{a,b)^i M aTbf' (2) 

tt' mm' 

where n a ^ is a normalization factor, FT is a func- 
tion of the power spectra Cf , which include both CMB 
lensing and primary anisotropy contributions, and ab = 
TT, TE, EE, EB, TB^tiie BB case is excluded because 
the method of Ref. [29| is only valid when the lens- 
ing contribution is negligible compared to the primary 
anisotropy, assumption that fails for the B modes in the 
case of Planck. 

It is possible to combine five quadratic estimators into a 
minimum variance estimator; the noise on the deflection 
field power spectrum Cf d produced by this estimator can 
be expressed as: 

N t d = V ( N aba'b'\-1 ■ ( 3 ) 

l^aa'bb' \ ly i ) 

A publicly available routine 

(http : //lesgourg. web. cern. ch/lesgourg/ codes .html) 

allows to compute the minimum variance lensing noise 
for the Planck experiment. At the same URL a full-sky 
exact likelihood routine is available and we use this in 
order to analyze our forecasted datasets, which include 
the lensing deflection power spectrum. 



We create a full mock CMB dataset with noise proj 
erties consistent with those expected for the Planck [1 
experiment (see Tab. [I]). 

Experiment Channel FWHM AT/T 



Planck 70 


14' 


4.7 


100 


10' 


2.5 


143 


7.1' 


2.2 


fsky = 0.85 







TABLE I. Planck-like experimental specifications. Channel 
frequency is given in GHz, the temperature sensitivity per 
pixel in fiK/K, and FWHM (Full- Width at Half-Maximum) 
in arc-minutes. The polarization sensitivity is assumed as 
AE/E = AB/B = V2AT/T. 



B. Galaxy weak lensing data 

We simulate future galaxy weak lensing data assum- 
ing the specifications for the Euclid weak lensing survey 
(see Table |TTJ) . This survey will observe about 30 galax- 
ies per square arcminute from redshift z = 0.5 to z = 2 
with an uncertainty of about a z = 0.05(1 + z) (see [l7j|). 
Using these specifications we produce mock datasets of 
convergence power spectra, again following the procedure 

of [13. 

The la uncertainty on the convergence power spec- 
trum {P{t}) can be expressed as (301 ] : 

y (2£+ l)f sky Ae \ n ga i J 
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n ga i(arcmin 2 ) redshift Sky Coverage y r 

(square degrees) 



30 0.5 < z < 2 15000 0.22 



TABLE II. Specifications for the Euclid like survey consid- 
ered in this paper. The table shows the number of galaxies 
per square arcminute (n ga i), redshift range, sky coverage and 
intrinsic ellipticity (7^ ms ) per component. 



where Ai is the width of the £-bm used to generate 
data, in our analysis we choose At = 1 for the range 
2 < £ < 100 and A e = 40 for 100 < £ < 1500. As 
at high £ the non-linear growth of structure is more 
relevant, the shape of the non- linear matter power 
spectra is more uncertain [31]; therefore, to exclude 
these scales, we choose £ m ax = 1500. We assume the 
galaxy distribution of Euclid survey to be of the form 
n(z) cx z 2 exp(— (z/zo) 1 ' 5 ) (see jl7j]). where z is set 
by the median redshift of the sources, zo — z m /lAl 
with z m = 0.9. Although this assumption is reasonable 
for the Euclid survey, the parameters that affect the 
shape of the distribution function may have strong 
degeneracies with some cosmological parameters as the 
matter density, as and the spectral index (32| . 



C. Analysis method 

We perform a MCMC analysis based on the publicly 
available package cosmomc [35[ with a convergence diag- 
nostic using the Gelman and Rubin statistics. 

The set of cosmological parameters that we sample is as 
follows: the baryon and cold dark matter densities ilt,h 2 
and ft c h 2 , the Hubble constant Ho, the scalar spectral 
index n s , the overall normalization of the spectrum A s 
at k = 0.002 Mpc -1 , the optical depth to reionization r, 
and, finally, the variation of the fine structure constant 
parameter a/ao . In our analysis we adopt flat priors on 
these parameters. 

We consider two cases. In a first run we assume a/ao — 
1 in the fiducial model and we investigate the constraints 
achievable on a and on the remaining parameters using 
the future simulated datasets. 

We then consider a fiducial model with a variation 
in a such that a/ao = 0.996, in principle detectable 
with the future data considered, and analyse the new 
dataset wrongly assuming a standard ACDM scenario 
with a/ao = 1. This analysis allow us to investigate how 
wrongly neglecting a possible variation in a could shift 
the best ?t cosmological parameters. 

In particular, since a variation in a essentially affects 
the recombination scenario at CMB decoupling, this ex- 
ercise is useful in understanding at what level of accuracy 
the recombination process should be computed in order 
to avoid a biased estimate of the main cosmological pa- 



rameters. 



III. RESULTS 

In Table IIIII we show the MCMC constraints on cos- 
mological parameters at 68% c.l. from our simulated 
dataset, obtained assuming a fiducial model with a/ao = 
1 

We consider two cases: a standard analysis where 
a/ao = 1 an d an analysis where also a/ao is varied. This 
is important in order to check at what level adding to 
the analysis one extra parameter affects the constraints. 
Moreover, in order to better quantify the improvement 
from the Euclid data we also report the constraints ob- 
tained using just the Planck data alone. 





Planck 


Planck+Euclid 


Model 


Varying a/ao 


a/ao = 1 


Varying a/ao 


a/ao = 1 


Parameter 










A(Q b h 2 ) 


0.00013 


0.00013 


0.00011 


0.00010 


A(Q c h 2 ) 


0.0012 


0.0010 


0.00076 


0.00061 


A(r) 


0.0043 


0.0042 


0.0041 


0.0029 


A(n 3 ) 


0.0062 


0.0031 


0.0038 


0.0027 


A(log[10 10 .4 a ]) 


0.019 


0.013 


0.0095 


0.0092 


A(H ) 


0.76 


0.43 


0.34 


0.31 


A(fi A ) 


0.0063 


0.0050 


0.0034 


0.0033 


A(a/a ) 


0.0018 




0.0008 





TABLE III. 68% c.l. errors on cosmological parameters from a 
first analysis made assuming a fiducial model with a/ao — 1. 



As we can see from the Table, the Euclid data im- 
proves the Planck constraint on a/ao by a factor ^2.6. 
This is a significant improvement since for example, a 
~ 2(7 detection by Planck for a variation of a could be 
confirmed by the inclusion of Euclid data at more than 5 
standard deviation. Moreover the precision achieved by a 
Planck+Euclid analysis is at the level of ~ 5 x 10~ 4 , that 
could be in principle further increased by the inclusion 
of complementary datasets. 

It is interesting to see what is the impact of a varia- 
tion in the fine structure constant in the estimate of the 
remaining cosmological parameters. There is a high level 
of correlation among a/ao and the parameters Hq and 
n s when only the Planck data is considered. This is also 
clearly shown in Figs. [1] and where we plot the 2-D like- 
lihood contours at 68% and 95% c.l. between a/ao, n s 
and Ho. Namely, a larger/lower value for a is more con- 
sistent with observations with a larger/lower value for Ho 
and a lower/larger value for n s . These results are fully 
in agreement with those reported in [? ] . 

When Planck and Euclid data are combined, the de- 
generacy with H is removed, yielding a better determi- 
nation of a. However the degeneracy with n s (see Figj2j) 
is only partially removed. This is mainly due to the fact 
that the n s parameter is degenerate with the reionization 
optical depth r, to which Euclid is insensitive. 
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0.1 - 




0.08 - 




0.994 0.996 0.998 



1 

a/a 



1.002 1.004 1.006 



0.075 L 



0.994 0.996 0.998 



1 

a/a 



1.002 1.004 1.006 



FIG. 1. 2-D constraints on a and Ho using Planck data (blue 
contours) and Planck+Euclid data (red contours). 



FIG. 3. 2-D constraints on a and r using Planck data (blue 
contours) and Planck+Euclid data (red contours). 



0.98 
0.975 

0.97 
0.965 

0.96 
0.955 

0.95 
0.945 




0.1 




0.994 0.996 0.998 



CC/CL 



1.002 1.004 



0.945 0.95 0.955 0.96 0.965 0.97 0.975 0.5 



FIG. 2. 2-D constraints on a and n 3 using Planck data (blue 
contours) and Planck+Euclid data (red contours) . 



FIG. 4. 2-D constraints on n a and r using Planck data (blue 
contours) and Planck+Euclid data (red contours). 



In fact, as can be seen in FigO using Euclid with 
Planck highlights a previously hidden degeneracy be- 
tween a/ao and t; both these parameters do not af- 
fect the convergence power spectrum, thus they are not 
measured by Euclid, but they are both correlated with 
other parameters, such as n s (Figf2] and FigfJJ, whose 
constraints are improved through the analysis of weak 
lcnsing. This improvement on n s allows to clarify the 
degeneracy between a/ao and r. From this discussion is 
clear that a better determination of the optical depth r, 
through, for example, future measurements of the CMB 
polarization, would further improve the constraints on a 
and other parameters. 

Clearly, as we can see from Table IIIII when a variation 
of the fine structure constant is considered in the analysis, 
the gain of including Euclid is significantly reduced in the 
constraints of n s and r. As we can see also from Table 
1, the errors on n s and r are increased by ~ 40% and 
~ 57% respectively when a is varying respect to the case 
when a is fixed to the standard value. 

This however should not bring us to the decision of not 



varying a in the future analyses of those datasets. What 
indeed could happen if the true value of a is different from 
the standard one and we perform an analysis wrongly 
assuming a/ao = 1 ? 

To answer to this question we have also analysed a 
mock dataset generated with a/ao — 0.996 (a varia- 
tion in principle detectable with this future data) but 
(wrongly) assuming a standard value (a/ao — 1)- 

The results, reported in Tab II VI show a consistent and 
significant bias in the recovered best fit value of the 
cosmological parameters due to the strong degeneracies 
among a/a and the Hubble constant H , the spectral 
index n s and the matter energy density Q m parameters. 

Note, from the results depicted in Figures [TJ [2] and 
Figures [5] and [5] and also from the results in Tab. II VI that 
the shift in the best fit values is, as expected, orthogonal 
to the direction of the degeneracy of a/ao with these 
parameters. For example, lowering a damps the CMB 
small scale anisotropics. As we can from[T]this effect can 
be compensated by increasing n s . Therefore a fiducial 
model with a lower value for a mimics a lower spectral 
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1 IdllL^iV 


-"RnrliH 


Hi Hnnal 


IA/n-1 








VtllllCS 




lVi.tJU.Cl. 


/ >. / — 1 

(JC I U.Q 1 


varying £; 






P 3x clyii c t. c r 












U.UZZOZ HI U.UUU1U 


n 099^0 -4- n nnm i 


U.UZZUO 


9 7 


n c h 2 


0.1129 ± 0.00059 


0.1106 ± 0.00078 


0.1109 


3.4 


T 


0.075 ± 0.0025 


0.088 ± 0.0041 


0.088 


5.3 


n s 


0.950 ±0.0028 


0.964 ± 0.0039 


0.963 


4.6 


Ho 


71.8 ±0.30 


71.0 ±0.33 


71.0 


2.7 


n A 


0.737 ±0.0032 


0.736 ± 0.0034 


0.735 


0.6 


0"8 


0.801 ± 0.0009 


0.803 ± 0.0010 


0.804 


3.3 



TABLE IV. Best fit values and 68% c.l. errors on cosmological 
parameters for the case in which a fiducial model with a/ao = 
0.996 is fitted wrongly neglecting a variation in a. The last 
column shows the absolute value of the difference between 
the best-fit value estimated fixing a/ao = 1 and the fiducial 
value, relative to the la error. 



index n Sl as shown in Tab. IIV1 

These results show that even for a small variation in 
a/ao ~ 0.5%, the best fit values recovered assuming that 
there is no variation in a can be at more than 95% c.l. 
away from the correct fiducial values, and may induce a 
significant underestimation of n s , r and as and an over- 
estimation of Hq. In the last column in Tab. II VI we show 
the difference between the wrong value estimated fixing 
a/ao — 1 an d the fiducial value, relative to the la error. 
We note that also other parameters, as fl c h 2 and ilfc/i 2 , 
have significant shifts. 

When a variation of a is considered, the correct fiducial 
values are recovered, however at the expenses of less tight 
constraints. 

Future analyses of high precision data from Euclid and 
Planck clearly need to consider possible deviations from 
the standard recombination scenario in order to avoid 
a significantly biased determination of the cosmological 
parameters. 



IV. CONCLUSIONS 

In this brief paper we have evaluated the ability of a 
combination of CMB and weak lensing measurements as 
those expected from the Planck and Euclid satellite ex- 
periments in constraining variations in the fine structure 
constant a. We have found that combining the data from 
those two experiments would provide a constraint on a 
of the order of Aa/a = 8 x 10~ 4 , significantly improving 
the constraints expected from Planck. These constraints 
can be reasonably futher improved by considering addi- 
tional datasets. In particular, accurate measurements of 
large angle CMB polarization that could provide a bet- 
ter determination of the reionization optical depth will 
certainly make the constraints on a more stringent. 

Moreover, we found that allowing in the analysis for 
variations in a has important impact in the determina- 
tion of parameters as n s , Hq and r from a Planck+Euclid 



0.975 




0.07 0.075 0.08 0.085 0.09 0.095 0.1 0.105 



FIG. 5. 2-D constraints on n s and r using a fiducial model 
with a/ao ~ 0.996, fitting it with a fixed a/ao = 1 (blue 
contours) and with a/ao aloowed to vary (red contours). 
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0.799 0.8 0.801 0.802 0.803 0.804 0.805 0.806 0.807 

FIG. 6. 2-D constraints on as and Ho using a fiducial model 
with a/ao = 0.996, fitting it with a fixed a/ao = 1 (blue 
contours) and with a/ao aloowed to vary (red contours). 



analysis. We have shown that a variation of a of about 
0.4% can significantly alter the conclusions on these pa- 
rameters. 

Changing the fine structure constant by 0.4% shifts 
the redshift at which the free electron fraction falls to 
x e = 0.5 by about ~ 1% from z* = 1275 to = 1262. 
An unknown physical process that delays recombination 
as, for example, dark matter annihilation (see e.g. }34|). 
may have a similar impact in cosmological parameter es- 
timation. Our conclusions can therefore be applied to the 
more general case of a modified recombination scenario. 
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